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We argue that the full conversion of a hadronic star into a quark or a hybrid star occurs within 
two different regimes separated by a critical value of the density of the hadronic phase rTh- The first 
stage, occurring for rih > nJi, is characterized by turbulent combustion and lasts typically a few ms. 

During this short time-scale neutrino cooling is basically inactive and the star heats up thanks to 
the heat released in the conversion. In the second stage, occurring for uh < nK, turbulence is not 
active anymore, and the conversion proceeds on a much longer time scale (of the order of tens of 
seconds), with a velocity regulated by the diffusion and the production of strange quarks. At the 
same time, neutrino cooling is also active. The interplay between the heating of the star due to the 
slow conversion of its outer layers (with densities smaller than TTh) and the neutrino cooling of the 
forming quark star leads to a quasi-plateau in the neutrino luminosity which, if observed, would 
possibly represent a unique signature for the existence of quark matter inside compact stars. We 
will discuss the phenomenological implications of this scenario in particular in connection with the 
time structure of long gamma-ray-bursts. 

PACS numbers: 21.65.Qr,26.60.Dd 


I. INTRODUCTION 

The Bodmer-Witten hypothesis on the absolute stability of strange quark matter [I|, stimulated many interesting 
investigations on the possible existence of compact stars entirely composed by this kind of matter Si or small 
nuggets of strange quark matter which would propagate in the Universe as cosmic rays S|. The exothermic process 
of conversion of ordinary nuclear matter into strange quark matter has been studied for the first time in Ref. Q 
where it has been modeled as a slow combustion by means of a one dimensional stationary reaction-diffusion-advection 
equation for the strange quarks concentration. This is a kinetic theory calculation in which the microphysical processes 
occurring within a finite width combustion zone are taken into account and which allows to determine the velocity of 
the conversion as a function of the quarks diffusion coefficient and the rate of conversion of down quarks into strange 
quarks. It turns out that typical values for the burning velocities are within 10^ — cm/sec for a quark chemical 
potential Hq ^ 300 MeV and a temperature of the quark phase T ^ 10 MeV. 

The limit of this kinetic theory approach is that it does not allow to take into account possible macroscopic collective 
flows and hydrodynamical instabilities driven by pressure and density gradients between the fuel and the ashes fluids. 
Gravity would also play an important role in the dynamics of the conversion. For a complete treatment of the 
problem one would have to couple the equations of hydrodynamics (i.e. the equations of conservation of baryon 
number, momentum and energy) and the equation of conservation of chemical species (which includes the diffusion 
and the reaction rates within the combustion zone) in multidimensional numerical simulations, see [3)0- The width of 
the combustion zone S can be estimated by the simple relation 6 ^ \/Dt Q where D is the quark diffusion coefficient 
{D ^ 10“^cm^/sec for pLq ^ 300 MeV and T ~ 10 MeV @) and r is the inverse of the rate of conversion of down 
quarks into strange quarks (r ^ 10“® sec for ^ 300 MeV 0)- One obtains d ^ 10 ® cm. Glearly, it would 
be numerically unfeasible to resolve such a small length scale within a numerical simulation aiming at studying a 
compact star whose size is of the order of ten km. 

A similar problem exists in the context of numerical simulations of type la Supernovae where length scales from 10“^^ 
to 10® cm characterize the physical system m and an alternative scheme has been devised: being the combustion 
zone much smaller than the size of the system one can assume that it is actually an infinitely thin layer and it can be 
treated as a surface of discontinuity which separates the ashes from the fuel. In this scheme, called flamelet regime [3^ . 
one has to impose the Hugoniot jump conditions to relate the thermodynamical variables of the fluid at both sides of 
the discontinuity. This approach, implemented in 3-l-lD, is one the most used in the context of type la Supernovae and 
it has demonstrated the crucial role played in such explosive events by the hy drodynamical instabilities, specifically 
the Rayleight-Taylor instability and the Landau-Darrieus instabilities [H, [l^ , which turns the laminar combustion 
into a much faster turbulent combustion. In this sense the flamelet approximation is very sensible: since the burning 
velocity is strongly enhanced by turbulence, the importance to know the exact value of the laminar velocity (governed 
by the microphysics of the combustion zone) is subordinate. 

The same method has been adopted for studying the conversion of nuclear matter into strange quark matter in 


2 


semi-analytical calculations [IM3- This process is most likely a deflagration (subsonic combustion) and not a 
detonation (supersonic combustion) although its velocity is substantially increased, up to two orders of magnitude, 
by the development of hydrodynamical instabilities. A common finding within this approach is that at some critical 
density this rapid combustion stops and the neutron star cannot fully convert into a strange quark star. These results 
have been recently confirmed in 3-l-lD hydrodynamics numerical simulations: the star is converted on a time scale of 
ms thanks to hydrodynamical instabilities but a sizable fraction of the star, of the order of few O.IMq, is left unburnt 
mill- A natural question that we will address in this paper arises: does the combustion really stop or does it 
proceed on a longer time scale? 

After clarifying the reason for which, in a purely hydrodynamical approach, the turbulent combustion stops (the 
critical density being given by the condition proposed by Coll in Ref. 0) , we will model the subsequent slow conversion 
process (which is driven by strangeness diffusion and weak interaction processes). We will demonstrate that in such 
a two steps scenario, a first stage of fast turbulent combustion and a second stage of slow laminar combustion, one 
could obtain a noticeable imprint on the neutrino signal emitted from the formation of a strange quark star. 

The paper is organized as follows: in Section II we will review the procedure adopted for treating the combustion 
within a purely hydrodynamical approach and we will explain the meaning of the Coil’s condition. In Section III we 
will model the slow combustion process and provide some example of the neutrino signals released by a strange star 
at birth. We draw our conclusion in Section IV. 


II. TURBULENT HYDRODYNAMICAL COMBUSTION 


We review here the way the combustion is treated if the combustion zone is so thin to be considered as a surface 
of discontinuity, usually called flame front. The generalization of the classical combustion theory within relativistic 
(magneto-)hydrodynamics has been presented by Coll and Anile in Refs.0, [2l|. Denoting with pi, Cj, Ui and 
Xi = (ei+pi)lnl the pressure, energy density, baryon density and dynamical volume of fluid i, the so-called “condition 
for exothermic combustion”, we will name it “Coil’s condition”, for the conversion of fluid 1 into fluid 2 reads: 
ei(p, X) > e 2 (p, X), i.e. at fixed pressure and dynamical volume, the energy density of fluid 1, the fuel, must be 
larger then the one of fluid 2, the ash, see also 22| for the case of classical hydrodynamics. Introducing the enthalpy 
density Wi = et +Pi, the Coil’s condition can be equivalently expressed as wi{p,X) > W 2 ip,X). As we will show, this 
condition is necessary in the case of a detonation, while in the case of a deflagration it determines the appearance 
of hydrodynamical instabilities dramatically accelerating the combustion process. To our knowledge this issue is not 
discussed in any standard textbook of hydrodynamics nor in any research paper. Our result also implies that Coil’s 
condition corresponds to the request of exothermicity in the case of detonations while its violation does not exclude 
the possibility of slow processes of combustion. 

Let us consider the two fluids to be hadronic matter and quark matter. For the hadronic matter we adopt a generic 
equation of state at zero temperature eh{ph,Xh) (this situation corresponds to the case of the conversion of a cold 
hadronic star) while for quark matter we consider the simple case of the equation of state of the MIT bag model with 


massless quarks as in Ref. 


The relation between energy density and pressure in this case reads: Cq = 3pq + 4:B 


where B is the usual bag constant (notice that in the case of massless quarks the energy density is a function of only 
p and not of X). The more general case of a polytrope is discussed in the appendix A, see also 0. 

Similarly to the case of the discontinuity associated with a shock wave, also in the case of the flame front, one has to 
impose the continuity equations for the fluxes of baryon number (or mass flux), momentum and energy. By indicating 
with j the number of baryons ignited per unit time and unit area of the flame front, one obtains the following relations 
between the thermodynamical quantities of the hadronic fluid and of the quark fluid: 


nhUh 

{pq-ph)/{Xh-Xq) 
'C^hiPhj XhfX^h Wqi^Pq, Xq)Xg 


UqUq = j 

f 

{Ph -Pq){Xh +Xq) 


( 1 ) 

( 2 ) 

( 3 ) 


the last equation corresponding to the so-called relativistic detonation adiabat. Uh and Uq are the four-velocities of 
hadronic and quark matter in the flame front rest frame. Starting from hadronic matter in the state A: ph = pa 
and Xh = Xa and for a given value of j, Eqs. 1-3 allow to determine the final state B of quark matter, pq = ps 
and Xq = Xb which lies on the detonation adiabat. In particular, the second equation represents a line in the (p,X) 
plane passing through A and with angular coefficient equal to —p. The intersections of this line with the detonation 
adiabat allow to find the state B of quark matter. The baryon number flux j, or equivalently the flame front velocity 
with respect to one of the two fluids Ui, in general cannot be expressed in terms of the thermodynamical variables 
of the states A and B. It is instead related to the specific microscopic properties of the chemical reactions of the 
combustion, the heat transfer and the diffusion of chemical species across the flame front and therefore it must be 
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FIG. 1: Illustrative plot of the detonation adiabat in the case in which the Coil’s condition is fulfilled (left panel) or not (right 
panel). A, B, B’ indicate respectively the initial hadronic state and two possible final states for the qnark phase. O and O’ are 
the Chapman-Jouget points. 


determined within a kinetic theory approach such as the one of Ref. @ . We remark however that hydrodynamical 
instabilities can substantially affect the flame front velocities as we will discuss later on. 

We will now show that Coil’s condition determines the position of the initial state A with respect to the detonation 
adiabat. To this purpose we will consider a point on the detonation adiabat having the same value of X of the initial 
state A {Xa = Xb) and we will compute the corresponding value of the pressure pb (this procedure corresponds to 
the limit j —>■ oo): in this way we will find if A lies below or above the detonation adiabat. Coil’s condition, written 
in terms of the enthalpy density, reads [l^: A(p, X) = eh{p,X) — eq{p, X) = Whip, X) — Wq{p, X) > 0 - can write 
Eq.3, adding and subtracting CqipA, Xa), as: 


Whip A, X a) - WqipB,XA) 
ehiPA,XA) +PA- eqipB,XA) -PB+ eqipA,XA) - eqipA,XA) 

Mpa, Xa) + CqipA, Xa) - tqipB, Xa) 
A(pa, Xa) + ipA + 45 - Sps - 45 
AipA,XA) = 2 (pb - pa) 


2pA - 2pB 
2 pA - ‘ 2 -Pb 
Pa - Pb 
PA -PB 


(4) 


Therefore if AipA,XA) > 0, i.e. if the Coil’s condition is fulfilled, then pB > PA which means that the initial state 
of the hadronic phase A lies in the region of the (p,X) plane below the detonation adiabat. In turn, this implies that 
there are two specific values of j, jo and jo', for which the lines passing through A are tangent to the detonation 
adiabat. The two points of tangency are the Chapman-Jouget points, see Fig.l. Following the standard treatment, 
one could obtain strong and weak detonations and strong and weak deflagrations depending on the specific values of 
j and on the boundary conditions of the problem. The point O corresponds to the Chapman-Jouget detonation and 
it is the only possible realization of detonation in a physical system, such a compact star, in which no external force 
is producing the shock wave, seeQ. Moreover a detonation taking place in a compact star can be assimilated to a 
detonation in a closed pipe and also in this case it can take place only at the Chapman-Jouget point Q. Clearly, if 
the Coil’s condition is not fulfilled, the state A lies in the region of the (p,X) plane above the detonation adiabat (see 
Fig.l right panel) and there are no Chapman-Jouget points in this case. Detonation is therefore excluded when the 
Coil’s condition is not fulfilled. 

Let us discuss now how the Coil’s condition is related to the deflagration regimes. We consider the simplest case 
of a slow combustion, i.e. a process in which the velocities Vh and Vq are much smaller than the sound velocities Ch 
and Cq of the two fluids (this case is particularly interesting for the combustion of hadronic stars into quark stars for 
which the laminar velocities, found in are much smaller than the sound velocities). By using Fqs.(l-3) one finds in 
this regime that Ph = Pa = Pq = Pb' and waI^a = wb'I^b' or equivalently (e^ +PA)l'nA = i^B' +PA)/nB', he. the 
enthalpy per baryon is conserved during the combustion (see (9| for the case of non-relativistic hydrodynamics). 
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FIG. 2: Difference between the energy density of the hadronic phase and the quark phase as a function of the baryon density. 
The solid black and red lines correspond to the LS180 and GM3 equations of state with only nucleons. The green lines 
correspond to the SFHo model with only nucleons (solid line), deltas (dotted line) and deltas and hyperons (dashed line). The 
Coil’s condition is fulfilled within the density window for which this difference is positive. 


One sees immediately that in this case the Coil’s condition implies, see Fig.l (left panel), that Xg > Xa i.e. 
(e_B' +pA)/n\, > (cA +PA)/nA which together with the conservation of the enthalpy per baryon implies ub' < nA- 
Moreover, from nA{eB' +Pa) = nB'{eA +Pa) < nA(fiA +Tt) on obtains bb' < ba- Thus the quark phase is produced 
with baryon density and energy density smaller than the one of the hadronic phase. This fact opens the possibility 
of obtaining the Rayleigh-Taylor instabilities. These instabilities arise, in the presence of gravity, if the gradient 
of the gravitational potential and the gradient of the energy density point in opposite directions (inverse density 
stratification). Indeed as shown in Refs. EH, the Rayleight-Taylor instabilities do occur during the conversion of 
an hadronic star and they substantially increase the efficiency of burning leading to time scales of the order of ms for 
the conversion of a big portion of the star. 

On the other hand, if the Coil’s condition is violated, the new phase is produced (again in the case of a slow 
combustion) with bb' > ba and therefore the burning can proceed but with velocities which are dominated by the 
diffusion and the rate of the chemical reactions and which are therefore much smaller than the velocities obtained 
during the turbulent regime. Moreover one can notice that since the new phase is produced with an energy density 
larger than the one of the fuel one cannot speak in this case of a deflagration (for which the inequality bb' < ba must 
hold true). 

One can separate the turbulent from the diffusive regime by finding the critical density of the hadronic phase, 
n/T, for which the Coil’s equality is satisfied Bh{p,X) = Bq{p,X). For this value of density the state A lies on the 
detonation adiabat. Moreover one also obtains that Ug = nT: thus pressure, energy density and baryon density are 
continuous across the interface of the flame front. The evolution of the flame for nh < Uh is not turbulent anymore. 

A simple way to understand the role played by the instabilities is to use the analytic model proposed in [1^ which 
allows to estimate the increase in the laminar velocities Ui due to turbulence in terms of the fractal dimension of 
the flame front. In particular, the fractal excess of the wrinkled front surface AD is proportional to the square of 
T = \ — Bb I Ba and it therefore vanishes when the Coil’s condition is met at ua = nT- The mean velocity Vrah of the 
front (with respect to the hadronic fluid) is significantly larger than the laminar velocity vih and it reads: 

‘Brah — (Aniax/'^min) (b) 

where Amax and Amin are two length scales which regulate the maximal and the minimal size of the wrinkle for which 
the velocity of the Rayleight-Taylor growing modes is larger than the laminar velocity, see As one can indeed 
notice in the simulations of [l8j| . the velocity of the flame front approaches the laminar velocity when the density is 
close to n^. Clearly a fraction of the star cannot be converted by means of a fast h ydr odynamical deflagration and 
actually a mass of a few O.IMq remains unburnt at the end of the turbulent regime |18l |. 
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FIG. 3: Enclosed gravitational mass and radius as a function of the baryon density for a l.SM© hadronic star before the 
turbulent conversion (black lines) and after the turbulent conversion (red lines). The black dashed line marks the appearance 
of hyperons: the seed of strange quark matter is formed at densities larger than this threshold. The red dashed line marks the 
density below which Coil’s condition is no more fulfilled and the turbulent combustion does not occur anymore. Below this 
density, the combustion proceeds via the slow diffusive regime. 


To provide some example of the results obtained by imposing the Coil’s condition, we display in Fig.2 the difference 
between the energy density of hadronic matter and of quark matter (again massless quarks for simplicity) for the 
following equations of state: the LS180 model from [2^, the GM3 model from [2^ and the SFHo model from [2^ 
also with the inclusion of delta resonances and hyperons [T7I| . Pure nucleonic equations of state provides two values 
of baryon density for which the Coil’s equation holds (as also found in [Hill). In general, the stiffer the hadronic 
equation of state the smaller the density window for which the turbulent hydrodynamical combustion can take place. 
One the other hand, when considering deltas and hyperons, within the SFHo model, the turbulent hydrodynamical 
combustion is always possible above the threshold for the formation of deltas (see dashed line). Notice that once a 
certain amount of hyperons is present in a hadronic star, the conversion is a necessary process in the scenario of the 
Witten’s hypothesis [^ . 


III. DIFFUSIVE REGIME 
A. Diffusion of strangeness and exothermicity 

As explained before, when the flame front reaches nt = nK the hydrodynamical instabilities responsible for the very 
fast conversion of the inner part of the star are not active anymore. The subsequent evolution of the system can be 
described by using the laminar burning velocity of the front vih computed in [y, : it corresponds to Vmh in the limit 

AD ^0 (see Eq. [g. Let us first review how the conversion process is described within microscopic kinetic theory 
approaches: we follow in particular the more recent treatment presented in M- The conversion is basically due to 
two simultaneous processes: diffusion of quarks within the combustion layer and flavor changing weak interactions 
among quarks (the process u + d ^ u + s being the most relevant for our case). One starts by introducing a position 
dependent strangeness unbalance a{x) = {nK{x) — n^)lnQ where the coordinate x = 0 defines the position of the 
conversion boundary, tik^x) = {nd{x) — ns{x))/2 and = nfc{x) for x —/ +oo (it would vanish for instance in the 
case of massless strange quarks). In the limit x —/ +oo one has bulk equilibrated strange quark matter with baryon 
density ng, while for a: < 0, a{x) = cost = a(0“) = un- If the hadronic phase is assumed to be made only of neutrons 
and if — 0 then qn = un/ nq with njv being the baryon density of neutron matter. In the following, for simplicity, 
we will adopt these conditions. Within the combustion zone, where both diffusion and chemical reactions are active, 
strange quark matter is out of beta equilibrium. There are two useful reference frames: the one in which the flame 
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front is at rest and the one in which quark matter is at rest. If vih and viq are the (laminar) velocities of hadronic 
matter and quark matter in the front frame, in the quark matter frame the velocity of the front is vif = —viq and 
v'lh = vih — viq. We remind that viq and vit are related by the ba^on flux continuity equation: viqUq = vihUh- 

By solving the steady state transport equation for a given in [y, llO| with the boundary conditions a(0''') = ag* = 
n^(0''')/nQ and a(x —>■ +oo) = 0 (notice that a{x) monotonically decreases from a: = 0 to a; —>■ +c») one obtains for 
the velocity of the front in the quark matter rest frame (in the suprathermal regime of Ref. 0): 


Vlf 


D 


a 


4 

Q* 


T 2aN{aN — OQ*) 


( 6 ) 


The quark diffusion coefficient D and the inverse rate of chemical reactions (for the process u + d ^ u + s) r are 
given by i,0: 


£1 = 0 . 


( Mg \2/3 

/ T \ 

V 300 MeV/ ' 

. 10 MeV ) 


-5/3 


cm^/sec, T = 1.3 X 10 


-9 


300 MeV 
Mg 


sec. 


(7) 


Notice that vif scales as ~ y-s/e. |.]^g more the star is heated up by the conversion process the slower the burning 
velocity. This fact will have important consequences for the thermal evolution of the star and the neutrino signal 
which will be discussed in the following. The value of the parameter og* is crucial: the larger ag* the larger vif. 
This quantity is just the boundary condition for the transport equation but its value cannot be taken arbitrarily large 
because one needs to impose the exothermicity of the process of combustion. In other terms, if the unbalance ag* is 
too large the front cannot move because there are not enough strange quarks to trigger deconfinement. Its maximum 
value ag“ is estimated in Appendix B. 

Let us now explain how we model the conversion during the diffusive regime. Once the initial state A of the hadronic 
phase is fixed, at a density Uh < n/T, one needs two equations to determine the state of the newly produced quark 
phase (for instance in terms of ^q and T). As in the case of a weak deflagration discussed before, since the velocities 
during the diffusive regime are small with respect to the sound velocities, we can set ph = Pa = Pq = Pb' (as also 
assumed in i0) and moreover Wh/rih = waIba = Wq/uq = wb'Ibb' [i3| Notice that these two equations match 
continuously a,t rih = Wh with the equations for the turbulent hydrodynamical regime. 

The process of conversion is exothermic if during the temporal evolution of the system some heat is released to 
the environment. Let us discuss how the heat per baryon generated by the conversion can be calculated. The 
equation of conservation of the enthalpy per baryon T^) = wb/bb^Pa^Tb) allows to find the value of 

Tb and if it turns out that Tb > Ta, the process is exothermic. By indicating with N the number of baryons 
composing the system, the total initial enthalpy is given (for uniform matter) by NwA/nA{pA,TA)- Because of 
cooling, asymptotically the system will reach again the initial temperature Ta and the total enthalpy of the system 
after the full conversion and after the cooling process will be NwB/nB{pA,TA)- Hence, the total heat Q released by 
the process of conversion and emitted into the environment (via neutrino cooling in the case of a compact star) is given 
by Q = N{wB/nBipA,TB) - wb/ nB{pA,TA)) = N{wA/nAipA,TA) - WB/nBipA,TA))- Thus, the difference between 
the enthalpy per baryon of the fuel and of the ashes calculated at the same pressure and temperature corresponds to 
the heat per baryon q released by the conversion: q = WA/nA{pA,TA) — wb/ nB{pA,TA)- 

A natural question concerns the point at which the conversion will stop. Let us indicate with phij") and Pq{r) the 
hadronic matter and quark matter pressure profiles inside the star (r is the radial coordinate). The radius of the 
star i?, is obtained by imposing p(r) = 0 whatever is the composition of the surface of the star. When the burning 
front is close to the surface of the star, the conservation of the enthalpy per baryon (for an initially cold star) implies 
eAlnA{TA = 0,_PA = 0) = eB/nB{TB > 0,pA = 0) > eB/nB{TB = 0,pA = 0). But eA/nA(TA = 0,PA = 0) = 
930 MeV (corresponding the energy per baryon of the iron nuclei composing the outer crust) and it is necessarily 
larger than eB/nB{TB = 0,pA = 0) under the hypothesis of absolutely stable strange quark matter. This means 
that the conversion is exothermic, and therefore it will continue, up to the surface of the star. Therefore, just by 
thermodynamics arguments, one would expect that the whole hadronic star converts into a quark star. On the other 
hand, the existence of a crystalline structure in the outer crust can significantly hinder the conversion process Q . It 
has been noticed in Ref. that the preheating of the crust due to the heat released by the conversion of the inner 
region leads to the dissociation of nuclei which greatly facilitates the burning into quark matter. In [2^ it has been 
estimated that the time needed to dissolve and convert the outer crust is of the order of few tens of ms. We will 
neglect in our calculation the conversion of the outer crust and consider only the, slower, conversion of the hadronic 
layer between and the neutron drip line density = 2.6 x 10“^fm“^. The conversion of the outer layer can be 
important e.g. when discussing short gamma-ray-bursts. We will consider this problem in a future paper. 
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B. Solving the diffusion eqnations 

We want now to discuss an example of the temporal evolution of the burning of a hadronic star during the diffusive 
regime. As a first ingredient, we need an initial profile for the star after the hrst stage of turbulent burning: this 
configuration is composed by hot quark matter for densities larger than n/j and by cold hadronic matter for densities 
smaller than rih. The equation of state of hot quark matter is computed as explained above by requiring that at fixed 
pressure, the enthalpy per baryon of the quark phase is equal to the one of the hadronic phase as in the case of a 
slow combustion. This assumption is not correct for the turbulent regime, during which the conversion velocity is 
significantly increased by the hydrodynamical instabilities; one can expect that the kinetic energy of the fluid flow is 
completely dissipated into heat once turbulence is over. 

We start from 1.5Mq hadronic star obtained with the SFHo model with the inclusion of deltas and hyperons. The 
central density of this stellar object is larger than the critical density for the formation of hyperons. In Fig.3, the black 
lines correspond to the enclosed gravitational mass (upper panel) and the radius (lower panel) as a function of the 
baryon density for the initial hadronic star configuration. The black dashed line indicates the onset of hyperons which 
are needed, in our scenario, to seed quark matter. The red lines correspond to the hybrid configuration, hot quark 
core and cold hadronic layer, after the turbulent regime. For the quark matter equation of state we have adopted the 
MIT bag model with massive strange quarks and with the inclusion of the perturbative QCD corrections (parameter 
set 1 of Ref. M)- Finally, the red dashed line located at n = Uh separates the turbulent from the diffusive regime. 
During the conversion the total baryonic mass, which turns out to be Mb = I.IIMq, is conserved. Notice that the 
total gravitational mass is also conserved during the turbulent conversion: this is due to the fact that cooling is not 
active during this stage and can be neglected. After the turbulent regime, a layer of hadronic matter of about 0.5Mq 
with a width of about 3km is left unburnt. This situation is exactly the one obtained also in the numerical simulation 
of Ref.[l^ and depicted in Figs 2 and 3 of that paper. 

Let us now introduce the two differential equations describing the propagation of the flame front and the thermal 
evolution of the star. We work in the reference frame of quark matter. Concerning the position of the flame front, by 
labeling with rf(t) its radial coordinate, one can write: 

=Vlf{fJ.q,T) ( 8 ) 

with the initial condition rf{0) = r where r is obtained from the baryon density profile by using the equation 
nh{r) = nK- 

For handling the thermal evolution of the star, in principle, one should couple Eql8] with a partial differential 
equation describing the heat transport through neutrinos and the heat source term related to the energy released by 
the conversion. Moreover one should also introduce the thermal conductivity which determines how the new release 
of energy is distributed within the star. The effect of gravity should also be considered to take into account the 
readjustment of the star during this stage. This is a very complicated task that we will not face in this work. We 
will adopt instead some simplifying assumptions that allow us to understand qualitatively how the conversion process 
proceeds and to obtain some order of magnitude estimates. 

First we can notice that the equations describing heat diffusion (and also the readjustment of the star) have been 
numerically solved in Ref. M- In that paper the diffusion regime was not discussed, but since strangeness diffusion 
is a rather slow process we can assume that the first few seconds of the thermal evolution of the star are dominated 
by the diffusion of the heat deposited during the rapid burning of its central re gion . We therefore assume that the 
temperature of the surface of the star and the neutrino luminosity evaluated in [l9| are the correct ones du ring the 
first seconds, even if the slow conversion of the outer layer is not taken into account. From Fig. 3 of Ref.fl^ one 
can also notice that after a few seconds the thermal profile inside the star flattens and that in particular the outer 
region (r > r) is almost isothermal after about 7-8 seconds. This first contribution to the neutrino luminosity can be 
approximated by using the simple formula 

L{t) = Qlre-*/^ (9) 

where Q is the total energy released during the rapid conversion and r is the time-scale of neutrino diffusion [30l| . 
The results of [l^ are well approximated by taking r ~ 3 s and Q ^ 8.5 x 10®^ erg. 

To simplify the description of the conversion of the outer layer we assume that: i) the conduction of heat is extremely 
fast (infinite thermal conductivity) in comparison with the burning velocity and therefore the heat generated by 
the conversion is distributed throughout the star (the very same assumption has been made in Q); ii) we assume 
that neutrinos are completely trapped and only the black body surface emission is considered for the cooling. The 
corresponding luminosity reads L = 21/8a{T/K)^4:TTr‘^ erg/s with Vg the radius of the neutrinosphere (we will 
assume that it is located at the interface between the inner crust and the outer crust where Uh = nd-) These two 
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FIG. 4: Neutrino luminosity associated with the burning during the diffusive regime of the combustion for three choices of 
the parameter The black line represents the luminosity obtained from the rapid combustion of the core as estimated by 

using EqEl 


assumptions imply that the temperature is uniform within the star and one can write a simple equation that expresses 
the energy conservation: 


C{T)— = -L{T) + Airr} j{rf, T) q{rf, T) (10) 

where C is the heat capacity of the star, L the neutrino luminosity and j = nhv'u^ is the number of baryons ignited 
per unit time and unit area. The thermodynamical variables nt, q and (which appears in are all functions of 
rf{t) and T{t). Concerning the heat capacity, we use C = 2 x 10^®M/Mq(T/10®) erg/K obtained in Ref. for a 
uniform density quark star or a hadronic star. 

By solving simultaneously Eqs. 15] and [TUI with initial conditions: r(0) = f, T(0) = To MeV (which is the temperature 
of the star for r > r after the turbulent regime and it is of the order of 5 MeV as found in [1^) we can calculate the 
time needed to complete the conversion of the star and the neutrino luminosity due to the conversion of the material 
left unburnt after the turbulent stage. In FiglU we show three cases corresponding to different values of {qn is 
fixed to one as in [^). We also display the curve of luminosity corresponding to Eq. 

As discussed above, we assume that the neutrino luminosity estimated in and approximated by Eq. [9] is close 
to the exact one during the first seconds. After some 10 seconds the luminosity estimated by taking into account the 
conversion of the outer layer becomes larger than the one obtained in [T^ where the outer layer remained unburnt. 
From that moment forward we assume that the estimate for the neutrino luminosity obtained by solving Eqs. [5] 
and [TU] is close to the exact one. In other words, we have separated the complicated problem of the heat diffusion 
and of the burning of the outer layer into two stages: the first one lasting a few seconds during which the neutrino 
luminosity is dominated by the heat dep osited during the rapid combustion of the core and is evaluated by solving 
the problem of heat transport as in [l9|| and a second stage, starting after a few seconds, during which the process 
of exothermic conversion of the outer layer provides enough heat to dominate the temperature of the star and the 
neutrino luminosity. During this second slow stage we have assumed the star to be isothermal. 

A quasi-plateau in the neutrino luminosity associated with the combustion of the hadronic layer is obtained (partic¬ 
ularly evident for the smallest value of This feature is a necessary consequence of the temperature dependence 

of the burning velocity: as the conversion proceeds, the temperature increases due to the release of energy and there¬ 
fore the velocity decreases. It is a self-regulating mechanism which rapidly leads to an almost constant velocity of 
burning and an almost constant luminosity of neutrinos. The process goes on until the whole star is converted. The 
kink appearing in the luminosity curves signals the end of the conversion: the following evolution is governed only by 
the cooling and the standard power law luminosity is obtained. 
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IV. CONCLUSIONS 


The conversion of a hadronic star into a strange quark star and the related neutrino emission can be divided into 
three different stages: 


• turbulent conversion of the inner part of the star, r < r, lasting a few ms. During this phase cooling is negligible 
due to neutrino trapping. 

• Diffusive conversion which burns the region r > r: 

a) during the first few seconds (at least for stars initially cold and having masses < I. 5 M 0 ) the neutrino 
luminosity is dominated by the cooling of the inner region r < r 

b) neutrino luminosity dominated by the cooling of the outer region r > r. This phase lasts a few tens of seconds 
and displays a quasi-plateau which originates from the burning velocity which is inversely proportional to the 
temperature. 


The quasi-plateau produced during the phase b) has a rather large luminosity ~ iQ5i-52 gj.gy'g fQj. gf gggonds 
and, if detected, would be a unique feature of the conversion of an hadronic star into a quark star. This very interesting 
possibility deserves a more refined treatment of the cooling process possibly via a neutrino transport code. 

The prolonged conversion of the star is very promising also in connection with the mechanism generating long 
gamma-ray-bursts within the protomagnetar model 31|. This model needs three crucial ingredients: high rotation 


frequency, high magnetic field and a significant neutrino emission (provided by the cooling of the protomagnetar) 
which, through ablation of the external layers of the compact star, allows to obtain the correct value of the Lorentz 
factor of the wind within which the gamma-ray-burst is generated. The neutrino emission caused by the conversion of 
an old hadronic star could lead to a long gamma-ray-burst if the other two requirements (i.e. high magnetic filed and 
high rotation frequency) are fulfilled. One possibility is given by a merger of a neutron star with a white dwarf. It has 
been shown that this process produces a spinning Thorne-Zytkow-like object with a low temperature, T ^ 10® K [^ . 
If large magnetic fields are generated, for instance via magnetorotational instabilities, the conditions for producing a 
gamma-ray-burst are fulfilled. The accretion of matter onto the hadronic star would trigger the conversion to a quark 
star and the expected neutrino luminosities are similar to the ones presented in Fig. 4. Such a gamma-ray-burst would 
be similar to a short gamma-ray-burst because it is associated with the merger of two compact stars but its duration 
would be comparable to the one of long gamma-ray-bursts. These features are in agreement with the analyses of 


GRB060614 [33l - l35j where it has been argued that this burst is not associated with a supernova. 


G.P. acknowledges financial support from the Italian Ministry of Research through the program “Rita Levi Montal- 
cini”. The authors gratefully acknowledge the COST Action MP1304 ’’NewCompStar” for supporting their networking 
and collaboration activities. 
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V. APPENDIX A 


We generalize the discussion presented in Sec.II, where we have adopted the quark matter equation of state with 
massless quarks, to the case in which the equation of state of fluid 2 is a generic polytrope given by: = auq+pq/ ( 7 — 1 ) 

([ 3 ^), Pq = kn'^, where 7 is the adiabatic index and 1 < 7 < 2 (the second inequality implying that the equation of 
state is causal at all densities [ 13 ). 

One can derive the following expression for the energy density as a function of Pq and Xq-. 


^q{Pq: Xq) 


Q!^(7 - 1) + 2pqXq + 0^7 “ 1 •\/a^ (T - 1) +^'lPqXq 

2^,(7-1) 


( 11 ) 


Let us fix Xb = X^ and assume X(pa,Xa) > 0. If pb > Pa^ the initial point A lies in the region of the (p,X) 
plane below the detonation adiabat. The detonation adiabat reads (adding and subtracting eq{pA, Xa))'- 


A{pa,Xa) = pa - pb + eq{pB,XA) - eq{pA,XA) ( 12 ) 

which after some manipulation and using Eg IIII reads: 

A{pa,Xa) = 2Aaa/7^ (\/«^(7- 1) +‘i'yPBXA - ^/a‘^ij-1) + AjpaXa^^ (13) 

+ (pb-Pa)- —7 (14) 

7-1 

Since 1 < 7 < 2 the sign of A{pa,Xa) clearly determines the sign of pb — Pa- Thus, if A(p^, A^) > 0, i.e. if the 
Coil’s condition holds true, the initial point A lies in the region of the (p,X) plane below the detonation adiabat. One 
sees again the Coil’s condition establishes the position of the initial state of the hadronic phase with respect to the 
detonation adiabat. [13 


VI. APPENDIX B 


In this section we aim at providing an estimate on the allowed values of the parameter ag* = (n.(j(0+) — 
ns(0'''))/(2nQ). Please notice that in this definition, appear quantities computed at the interface and a quantity 
computed for x —>■ + 00 . In particular we estimate which is the maximum value of the strangeness unbalance at 
the interface between hadronic matter and quark matter in order for the process of conversion to still proceed because 
of its exothermicity, see Fig.4 of Ref. We assume, as in Ref. [13 , that the deconfinement always takes place with 
maximal velocity i.e. with ag* = ag“. We consider the case of a bag-model equation of state with massless quarks 
and we limit our discussion to the case of vanishing pressure p. This case is particular relevant for the burnin g o f the 
most external layers of the hadronic star where the pressure approaches zero. Following the scheme of Ref. [38|, we 
introduce the parameter which is the ratio between the baryonic density of strange quarks and the baryon density: 
Xg = ns/3n. The Fermi momenta of up, down and strange quarks read: j^{n^rs) = kF{n,rs) = (1.5 7 r^n(l —rg))^/^, 

kFg{n,rg) = ( 37 r^(rsn))^/^ where we have assumed that up and down quarks have the same Fermi momenta (as in 
beta stable quark matter for the case of massless quarks). The energy density reads: 


e = 


Qkp 


34 ^ 

47r^ 


B. 


(15) 


The bag constant B varies in the range Bmin < B < Bmax where for B = Bmax the energy per baryon of strange 
quark matter is 930 MeV and for B = Bmin the energy per baryon of two flavor quark matter is 930 MeV. 

We hrst compute nq which corresponds to finding the baryon density n = nq for which beta stable quark matter 
(i.e. Xs = 1/3) has vanishing pressure. At this value of density the energy per baryon (e/n)/? is strictly smaller than 
930 MeV because we are assuming the hypothesis of absolutely stable strange quark matter. 

We now calculate the maximum value of 71^(0'*’) — 718 ( 0 “*"). We define n and fj as the values of xi and Xg for which 
the pressure p = 0 and the energy per baryon of non beta-stable quark matter (e/n)non -/3 = 930 MeV. rj represents 
therefore the minimum amount of strangeness for which non-beta stable strange quark matter has an energy per 
baryon equal to the one of iron. A quark phase with Cg > 77 would be favored with respect to nuclear matter and the 
conversion process would therefore be exothermic. 

^max jg related to fj and n by the following equation: 


^max 

O-Q* 


kpin, Xg) - kl,^{n,Xg) 
2n‘^nq 


( 16 ) 












11 


For each value oi B, r^, ri and uq are completely determined and thus also is fixed. 

We display the results in Fig. [S] In general < 1, its maximum value equals 1 only if one considers for x —>■ +oo 
non beta-stable 2 flavor quark matter with Ud = 2n„ (see i). We consider instead, as in beta stable strange 
quark matter for x —>■ -l-oo. One can show analytically that in this case the maximum value of is 1 . 53 / 4/2 - 0.667 
and it is obtained for B = Bmin [13|- 

The minimum value of Oq/® is zero and it is obtained a.t B = Bmax- Notice that for this value of B, - 00 . 

This basically implies that the parameter space (i.e. the range of B) has little room for the minimum value of 
In other words, it is very unlikely that ^ 0. In turns, this means that while the value of the velocity of conversion 
is quite uncertain it cannot be vanishingly small. 
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